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The Tsallis distribution has been used recently to fit the transverse momentum distributions of 
identified particles by the STAR collaboration [l|] at the Relativistic Heavy Ion Collider and by 
the ALICE 0] and CMS 0] collaborations at the Large Hadron Collider. Theoretical issues are 
clarified concerning the thermodynamic consistency of the Tsallis distribution in the particular case 
of relativistic high energy quantum distributions. An improved form is proposed for describing the 
transverse momentum distribution and fits are presented together with estimates of the parameter 
q and the temperature T. 

Keywords: Tsallis, Thermodynamics, Consistency, Heavy Ions 



I. INTRODUCTION 

The Tsallis distribution has gained prominence re- 
cently in high energy physics with very high quality fits 
of the transverse momentum distributions made by the 
STAR collaboration [l| at the Relativistic Heavy Ion Col- 
lider and by the ALICE 2\ and CMS [3] collaborations 
at the Large Hadron Collider. 

In the literature there exists more than one version of 
the Tsallis distribution [4, 5] and we would like to inves- 
tigate in this paper one version which we consider suited 
for describing results in high energy particle physics. Our 
main guiding criterium will be thermodynamic consis- 
tency which has not always been implemented correctly 
(see e.g. The explicit form which we use for the 

transverse momentum distribution in relativistic heavy 
ion collisions is: 

dN 
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PtTUt coshy 
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where px and itlt are the transverse momentum and mass 
respectively, y is the rapidity, T and /i are the tempera- 
ture and the chemical potential, the other variables are 
defined below. In the limit where the parameter q goes to 
1 this reproduces the standard Boltzmann distribution: 
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In order to distinguish Eq. ([TJ from the form used by 
the ALICE and CMS collaborations @, HI we will refer 
to Eq. ([1]) as the Tsallis-B parameterization. Note in 
particular the extra power of q on the right hand side. 
The motivation for preferring this form is presented in 
detail in the rest of this paper, see in particular section 
VII. 



Thermal models have been successful in describing par- 
ticle yields at different beam energies [9MTT| , especially in 
heavy ion collisions. These models assume the formation 
of a system which is in thermal and chemical equilib- 
rium in the hadronic phase and are characterized by a 
set of thermodynamic variables for the hadronic phase, 
most important among these are the chemical freeze-out 
temperature and baryon chemical potential. The decon- 
fined period of the time evolution dominated by quarks 
and gluons remains hidden: full equilibration generally 
washes out and destroys large amounts of information 
about the early deconfined phase. 



While the description of integrated particle yields is 
reasonably successful, more detailed descriptions, espe- 
cially of the transverse and longitudinal momentum dis- 
tributions call for additional dynamics. The transverse 
momentum distribution is often described by a combina- 
tion of transverse flow and a thermodynamical statistical 
distribution. With the Tsallis distribution this superposi- 
tion is not used and a very good fit can be obtained using 
the additional parameter q which describes the deviation 
from a Boltzmann distribution. In the limit where q — > 1 
one recovers the standard statistical Boltzmann distribu- 
tion. Whether this is ultimately the correct description 
or not remains to be seen. This paper is a contribution 
to the understanding of the use of the Tsallis distribu- 
tion in high energy collisions. It is not meant as giving 
a final answer to the correct dynamical theory of heavy 
ion collisions. 

In the next section we review the derivation of the Tsallis 
distribution by emphasizing the quantum statistical form 
and the thermodynamic consistency. 



2 



II. TSALLIS DISTRIBUTION FOR PARTICLE 
MULTIPLICITIES. 

Several generalizations of the standard Fermi-Dirac 
distribution 



FD / — 



1 + cxp ( 2=£ 



(3) 



to a Tsallis form have been proposed in the literature, 
some of these have been shown not to be thermodynami- 
cally consistent. In the following we use the Tsallis form 
of Fermi-Dirac distribution proposed in [H, Il2l - fl5j which 
uses 



FD 



(E) 



1 



1 + cxp, 



(4) 



where the function exp (x) is defined as 



exp g (x) 



[\ + {q-l)x] l/{q ~ l) if 
[l + (l-g)x] 1/(1 ~ 9) if 



x > 
x < 



(5) 



and, in the limit where q 
exponential: 



1 reduces to the standard 



lim exp (x) — > exp(a 

q— >1 H 



The form given in Eqs. ((H) and (|5|) will be referred to as 
the Tsallis-FD distribution. The Bose-Einstein version 
(given below) will be referred to as the Tsallis-BE distri- 
bution (l6j while the Boltzmann approximation will be 
referred to as Tsallis-B distribution. It should be noted 
that variations of the above have been presented previ- 
ously in the literature. These will not be considered in 
this paper. 

As is well-known, all forms of the Tsallis distribution 
introduce a new parameter q. In practice this parameter 
is always close to 1, e.g. in the results obtained by the 
ALICE and CMS collaborations typical values for the 
parameter q can be obtained from fits to the transverse 
momentum distribution for identified charged particles 
and are close to the value 1.1 (see below). The value 
of q should thus be considered as never being far from 
1, deviating from it by 20% at most. An analysis of 
the composition of final state particles leads to a similar 
result [13| for the parameter q. 

In the limit where q — > 1 the two forms coincide. Nu- 
merically the difference is small, as shown in Fig. (TTJ) for 
a value of q = 1.1. 

The Boltzmann approximation leads to the result [1, [B[ 



i + (q - i 



T 



(0) 



Note that we do not use the normalized ^-probabilities 
which have been proposed in Ref. [5| since we use here 



Fermi - Dirac q = 

Tsallis-FD 




FIG. 1. Comparison between the Fermi-Dirac and Tsallis-FD 
distributions as a function of the energy E, keeping the Tsallis 
parameter q fixed, for various values of the temperature T. 
The chemical potential is kept equal to one in all curves, the 
units are arbitrary. 



mean occupation numbers which do not need to be nor- 
malized. In the limit where q — > 1 all distributions coin- 
cide with the standard statistical distributions: 



lim n%(E) = n (E), 



lim n^ D {E) 



FD 



lim n% E {E) = n BE (E). 



(7) 
(8) 
(9) 



A derivation of the Tsallis distribution, based on the 
Boltzmann equation, has been given in Ref. (l8j . A com- 
parison between the and n FD distributions is shown 
in Fig. (JTJ) . For the Boltzmann approximation, the Tsal- 
lis distribution is always larger than the Boltzmann one 
if q > 1. Taking into account the large px results for 
particle production we will only consider this possibil- 
ity in this paper. As a consequence, in order to keep 
the particle yields the same, the Tsallis distribution al- 
ways leads to smaller values of the freeze-out temperature 
for the same set of particle yields [l7|. The Tsallis dis- 
tribution for quantum statistics has been considered in 
Ref. 



ition tor quantum 
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III. THERMODYNAMIC CONSISTENCY 

The first and second laws of thermodynamics lead to 
the following two differential relations [21| 



de = Tds -\ 


- udn, 


(10) 


dP = sdT - 


- ndfj,. 


(11) 


where e = E/V, s — S/V and n — N/V. Since these are 
total differentials, thermodynamic consistency requires 
that the following relations be satisfied 
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The pressure, energy density and entropy density are 
all given by corresponding integrals over Tsallis distri- 
butions and the derivatives have to reproduce the cor- 
responding physical quantities. For completeness, in 
the next section, we derive Tsallis thermodynamics us- 
ing the maximal entropy principle and discuss quantum 
g-statistics in particular Bose-Einstein and Fermi-Dirac 
distribution by maximizing the entropy of the system for 
quantum distributions. This follows partly the derivation 
of Ref . Q . We will show that the consistency conditions 
given above are indeed obeyed by the Tsallis-FD distri- 
bution. 



IV. QUANTUM STATISTICS 

The entropy in standard statistical mechanics for 
fermions is given in the large volume limit by: 



S 



FD 



FD 



+ (l-n FD )ln(l-n FD )] 



(16) 



where g is the degeneracy factor and V the volume of 
the system. For simplicity Eq. (fTB) refers to one particle 
species but can be easily generalized to many by sum- 
ming over all of them. In the limit where momenta are 
quantized, which is given by: 

s FD = - g £ N ln n » + (i - n «) Mi - n ^ > (1 7 ) 



For convenience we will work with the discrete form in 
the rest of this section. The large volume limit can be 
recovered with the standard replacement: 



d 3 p 



(18) 



Thegeneralization, using the Tsallis prescription, leads 

to dMl 

s t d = K N th + (1 - in 9 (i - m)] , (19) 



where use has been made of the function 

x x - q - 1 



m g (2;) 



(20) 



often referred to as q-logarithm. It can be easily shown 
that in the limit where the Tsallis parameter q tends to 
1 one has: 



lim ln q (x) = ln(x). 



(21) 



The maximization of the entropy p^|) will give the rij's 
their Tsallis-type form. If we use the explicit form of the 
q-logarithms we obtain 



(1 - m) - (1 - nrf 



1 



1 



c fd \ " 

i L \ * - / \ 1 - 

(22) 

In a similar vein, the generalized form of the entropy for 
bosons is given by 

s t e = -9^2 k ln « n - - c 1 + n ^ q ln «(i + "*)] ' ( 23 ) 



by using a similar method, we can express equation (|23j) 
as 



{i + ni )-{i + ni y 

9-1 



(24) 

In the limit q — > 1 equations (fT9|) and (f23|) reduce to the 
standard Fermi-Dirac and Bose-Einstein distributions. 
Further, as we shall presently explain, the formulation 
of a variational principle in terms of equation (|22[) allows 
to prove the general relation of thermodynamics. One of 
the relevant constraints is given by the average number 
of particles, 



(25) 



Notice the unusual power of q on the left-hand side. As 
it turns out, it is necessary to have this power of q since 
otherwise there is no thermodynamic consistency. 
Likewise, the energy of the system gives a constraint, 



(26) 



again, it is necessary to have the power q on the left-hand 
side as no thermodynamic consistency would be achieved 
without it. The maximization of the entropic measure 



4 



equation ([22]) under the constraints equation (|25|) and 
(f2"6")) leads to the variational problem. 



_S_ 

8m 



S FD + a{N _ £ n?) + ^ _ £ ^ 



= 0. 



(27) 

where a and /? are Lagrange multipliers associated, re- 
spectively, with the total number of particles and the 
total energy. Differentiating each expression in equation 
03 ' 



8m 



<?-i 



1 - n, 



9-1 



- 1 



(28) 



which leads to 

"dm ^ >~^1 



1 + m 



9-1 



,9-1 



(35) 



and by using equations (l3"5l) . (|2T)|) and ([3U|) in (jM)) . one 

gets 



1 + n., 



9-1' 



- ^ - a j = 0. 

(36) 

By rearranging equation (I36[) . one obtains the expression 
for the Tsallis-BE distribution [H, 
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and 
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(29) 



(30) 



then by substituting equation (|28p. (j2"9")l and ([3"U| into 
(l27l). we obtain 



qn 



9-1 



-1 



i - m 



which can be rewritten as 



9-1 
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9-1 



-1 



1 — n,- 
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- PEi - aj = 0. 

(31) 

= + a, (32) 



and, by rearranging equation J32J), we get 
^— ^ = [l + (g-l)C8Bi + a)]^ 



finally, we get the solution of generalized form of Fermi- 
Dirac distribution like this 

1 



[l + {q-\){pE l + a)] — + 1 
1 



[exp q (a + PEi)] + r 



(33) 



Which is the expression for the Tsallis-FD distribution 
referred to earlier in this paper 

Using a similar approach one can also determine the 
Tsallis-BE distribution [l6[ . Starting from the extremum 
of the entropy subject to two conditions one has: 



8n, 



= 0, 
(34) 



[l + {q-\){pE l + a)]— - I 
1 

[exp 9 ((£i-/i)/T)] -1' 



(37) 



where the usual identifications a = —\ijT and /3 = 1/T 
have been made. 



V. PROOF OF THERMODYNAMIC AL 
CONSISTENCY 

In order to use the above expressions it has to be shown 
that they satisfy the thermodynamic consistency condi- 
tions. To show this in detail we use the first law of ther- 
modynamics [2~fl | 



P 



-E + TS + (iN 
V ' 



(38) 



and take the partial derivative with respect to \x in order 
to check for thermodynamic consistency, it leads to 



dP 
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1 

V 



be as dN 
-ir + T ir + N + ^ir- 

o/i on op,, 

ET ( , , Ej — u\ dn q 
-— jfl + fe-l) 



T 



| Tgji - my- 1 dm 

q — 1 ctyt 
then, by explicit calculation 

dn\ _ qnj +1 
~d/7 ~ T 



(39) 



1 + («-!) 



Ei-n 



T 



dm _ nf 



1 + («-!) 



T 



-1+T^T 



and 



(1-ni) 



9-1 9-1 

= n; 



(g-l)(^-/x) 
T 



5 



Introducing this into equation (|39|) . yields 

dP 
du 



(40) 



which proves the thermodynamical consistency (|14[) . 
We also calculate explicitly the relation in equation 
can be rewritten as 



dE 

as 



n QT ' 9/i " 
_ 9T g/j dT 



I dS dp ' 
dT ' 9/i dT 



(41) 



since n is kept fixed one has the additional constraint 



leading to 



dT 
a/i 



Now, we rewrite (|4Tj) and (|42|) in terms of the following 
expressions 



^=E<^ ot' 



9/i 



^q-idrij 
dT' 



E<^ 



q-ldUj 



dS_ 

df 



E<? 



n?" 1 + (1 ~ n^" 1 
9-1 



as 



= E* 



9-1 



9-1 



drii 
'dT 1 



drij 
d\x ' 



and 



dn _ 1 

~ V 



dn _ 1 
djl ~ V 



E 



~dT 



By introducing the above relations into equation (1411) . 



the numerator of equation (|4Tj) becomes 
dE dE d/i m—^ q _ 1 dn i 
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Where the abbreviation 

Cij = {n i n ) q ^ 1 



drii drij dnj drii 
~dT~dn~ ~ dT dn 



(44) 



has been introduced. One can rewrite the denominator 
part of equation (|41l) as 



(42) 95 dS_djJ_ 
dT d[i dT 



-EH^ + a-^r 1 



. i 9nj 
djj, 
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where 



rip 1 + (1 - tn)'- 1 (^- m ) 
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T 



hence, by substituting equation (|4"3"]l and (|4"5")) in to (|41l) . 
we find 



as 

95 



E E i c k 



T- 



1,3 



(46) 



since Cy = 0, this finally leads to the desired result 

(47) 



'■j 



dE 
~dS 



= T. 



Hence thermodynamic consistency is satisfied. 

VI. BOLTZMANN APPROXIMATION 

Due to its practical relevance and importance we de- 
vote a section to the Tsallis-B distribution. In this case 
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the entropy is obtained from equation (|15l) by assuming 
the rii <C 1, this leads to 



w , q . 



(48) 



i=l 



The n$ are given explicitly as 



! + (<?-!) 



T 



(49) 



where rij denotes the number of particles in the ith en- 
ergy level with energy £7j. The maximum of the above 
entropy is looked for under the constraints imposed by 
fixing the total number of particles N and the total en- 
ergy of the system E, as given in equation (f^5|) and (f2l>| . 
As in the previous section, it should satisfy thermody- 
namic consistency which is given in equation (|12[) . The 
derivative of pressure w.r.t. becomes 



dP 
9/i 



1 
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dE 
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now, by using 



and 



dn\ a-idrii 



drii _ n\ 
d/7 ~ T~' 



By the above relations in equation ([5H)l , we recover equa- 
tion (HI). 



We now calculate the expressions needed in equations 
ID and dH) in terms of 
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then the numerator part of equation (|4Tj) become 



d£ as dfi 
df + ~djldf 



-\drij 
dT 



E^j 



_i drij 



'^2qE l (n i n J ) q 1 dj 



1,3 



3 



(51) 



Similarly, the denominator part of equation (|41[) can be 
written as 

.9-1" 



dS dS dfj, 



dT <9/z dT 
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where 
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I - qn 
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by combining the expressions in equation (|5ip and (|52[) 
into (HU, we find as before 

dE 



OS 



T. 
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It has thus been shown that the definitions of tempera- 
ture and pressure within the Tsallis formalism for non- 
extensive thermostatistics lead to expressions which sat- 
isfy consistency with the first law of thermodynamics. 



^ - V 
~d^~^ 



1 - qn 
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dn t 
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while the other partial derivatives are the same as previ- 
ously, by plugging the above relations into equation (1411) , 



VII. THERMAL FIT DETAILS 

The total number of particles is given by the integral 
version of (1251), 



N = gV 



d 3 p 
(2^)3 



1 + (<?-!) 



E- /2 
T 



9/(1-9) 



(54) 
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The extra power of q is necessary for thermodynamic 
consistency. The corresponding (invariant) momentum 
distribution is given by 



E §-p =9VE J^f 



1 + (q - 1 



E-fi 



1-2/(1-9) 



T 



(55) 



which, in terms of the rapidity and transverse mass vari- 
ables, becomes 



values quoted by the CMS collaboration Q . The analytic 
expression used in Refs. [2|, |3[ corresponds to identifying 



(59) 



an additional factor of the transverse mass on the right- 
hand side and a shift in the mass. 



VIII. CONCLUSIONS 
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At mid-rapidity y = and for zero chemical potential 
this reduces to the following expression 
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or, introducing the transverse momentum: 
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Fits using the above expressions based on the Tsallis-B 
distribution to experimental measurements published by 
the CMS collaboration Q are shown in Figs. 2, 3 and 4 
and are comparable with those shown by the CMS col- 
laboration but the resulting parameters are considerably 
different and are collected in Table I. The most striking 
feature is that the values of the parameter q are fairly 
stable around the value q w 1.11 for all cases consid- 
ered, whether 7 Tev or 0.9 TeV. The same cannot be 
said about the temperature T which is around 100 MeV 
with considerable deviations, it is however well below the 



In this paper we have presented a detailed derivation 
of the quantum form of the Tsallis distribution and con- 
sidered in detail the thermodynamic consistency of the 
resulting distribution. It was emphasized that an addi- 
tional power of q is needed to achieve consistency with the 
laws of thermodynamics Q. The resulting distribution, 
called Tsallis-B, was compared with recent measurements 
from the CMS collaboration Q and good agreement was 
obtained. The resulting parameter q which is a measure 



Particle 


T (MeV) 


q 


K% (0.9 TeV) 


92 


1.13 


K% (7 TeV) 


105 


1.15 


A (0.9 TeV) 


70 


1.11 


A (7 TeV) 


117 


1.12 


3" (0.9 TeV) 


44 


1.11 


ET (7 TeV) 


126 


1.11 



TABLE I. Fitted values of the T and q parameters for strange 
particles measured by the CMS collaboration 0] using the 
Tsallis-B form for the momentum distribution. The normal- 
ization has been adjusted. 



for the deviation from a standard Boltzmann distribution 
was found to be around 1.11. The resulting values of the 
temperature show a wider spread around 100 MeV. 

Whether or not the Tsallis distribution provides a valid 
interpretation of high energy collision data will need fur- 
ther theoretical work. 



[1] B. I. Abele v et al. ( STAR), Phys . Rev.| C75, 064901 

(2007), |arXiv:nucl-ex/0607033| 
[2] K. Aamodt et al. (ALICE Collaboration)(2011), 



[7] 

[8] 



|arXiv:1101.4110 [hejp^ex] 
[3] V. Khachatryan et al. ( CMS), IJHEPI 05, 064 (2011), 

|arXiv: 1102 4282 [hep-ex]) 
[4] C. Tsallis, |J.Statist.Physl 52, 479 (1988). 
[5] C. Tsallis, R. S. Mendes, and A. R. Plastino, Physica 

A261, 534 (1998). 
[6] F. Pereira, R. Silva, and J. Alcaniz, 



015201 (2007), arXiv:070 5.0300 [n ucl-th 
F. Pereira, R. Silva, and J. Alcaniz, P 
4214 (2009), |axXiv:0906.2422 [nucl-th] | 
J. M. Comoy; B. G~ Miller, 



'hys.Rev. C76, 



iys.Lett-1 A373, 



and A. R. 



Plasti no, |Phys.Lett.| A374, 4581 (2010), 

arXiv: 1006.3963 [cond-mat.stat-mech] 
[9] J. Cleymans, H. Oeschler, K. Redlich, and S. Wheaton, 

|Phys. Rev.| C73, 034905 (2006), [ar3SV:hep-ph/05iTu94l 
[10] A. Andronic, P. Braun-Munzinger, and J. Stachel, 

|Nucl. Phys.] A772, 167 (2006), [arXiv:nucl-th/0511071| 
[11] F. Becattini, J. Manninen, and M. Gazdzicki, Phys. Rev. 

C73, 044905 (2006), [arX"iv:0511092 [hep-ph]| 
[12] F. Buyukkilic and D. Demirhan, Phys.Lett. A181, 24 

(1993). 

[13] F. Pennini, A. Plastino, and A. R. Plastino, 

|Physics Letters A| 208, 309 (1995). 
[14] A. M. Tew eldeberhan, A. R. Plastino, and H. G. Miller, 

|Phys.Lett.| A343, 71 (2004). 



8 



[15] J. M. Conroy and H. Miller, |Phys.Rev.| D78, 054010 

(2008), |arXiv:0801.0360 [hep-ph]J " 
[16] J. Chen, Z. Zha ng, G. Su, L. Chen, and Y. Shu, 

|Physics Letters A 300, 65 (2002). 
[17] J. Cleymans, G. Hamar, P. Levai, and S. Wheaton, 

|J.Phys.| G36, 064018 (2009 ) |arXiv:0812.1471 [hep-ph]| 
[18] T. Bir o and G. Purcsel, |Phys.Rev.Lett.| 95, 162302 

(2005), |arXiv:hep-ph/0503204 [hep-ph]| 
[19] A. Teweldeberhan, H. Miller, and 

R. Tegen, |Int.J.Mod.Phys. E12, 395 (2003), 

|arXiv:0210011 [hep-ph]| 
[20] A. Plastino, A. Pl astino, H. Miller, and H. Uys, 

|Astrophys.Space Sci.| 290, 275 (2004). 
[21] S. R. de Groot, W. A. van Leeuwen, and C. G. van Weert, 

Relativtsttc Kinetic Theory (North Holland, 1980). 



— Tsallis-B fit 




I i I i I i I i L 

2 4 6 8 

P T (GeV) 



FIG. 2. Comparison between the measured transverse mo- 
mentum distribution for Kg as measured by the CMS collab- 
oration 3] and the Tsallis-B distribution as given by Eq. (|57[) 
using the parameters listed in Table I. 
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FIG. 3. Comparison between the measured transverse mo- 
mentum distribution for A as measured by the CMS collabo- 
ration [3[ and the Tsallis-B distribution as given by Eq. (|57[) 
using the parameters listed in Table I. 
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FIG. 4. Comparison between the measured transverse mo- 
mentum distribution for A as measured by the CMS collabo- 
ration |3| and the Tsallis-B distribution as given by Eq. (|57[) 
using the parameters listed in Table I. 



